Channel measurement method for large-scale antennas, and user terminal

ABSTRACT

Provided are a channel measurement method for large-scale antennas, and a user terminal. The method comprises: antennas at a base station side are divided into more than two groups in advance, a user terminal receives pilot signals transmitted from the antennas at the base station side, where the antennas in one group transmit the pilot signals in a frequency-division mode, and the antennas in different groups transmit the pilot signals in a time-division mode; and the user terminal performs a spatial channel estimation based on the pilot signals received from partial groups of the antennas, to obtain a channel measurement result for all the antennas.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of International Application No. PCT/CN2014/082734, filed on Jul. 22, 2014, which claims priority to Chinese Patent Application No. 201310312322.5, filed on Jul. 23, 2013, both of which are hereby incorporated by reference in their entireties.

FIELD OF THE TECHNOLOGY

The present application relates to the field of wireless communication technology, particularly to a channel measurement method for large-scale antennas and a user terminal thereof.

BACKGROUND

The main object of the next generation cellular communication system is to improve the system throughput, and MIMO (Multiple-Input-Multiple-Output, multiple-input-multiple-output) is an important technical solution to achieve the object in modern mobile communications. The technology of Large-scale antennas is further development of the MIMO technique, which improves the system throughput by deploying massive antennas at the base station side, as shown in FIG. 1.

In an MIMO system, the base station is required to transmit pilot signals to the whole cell at intervals, a UE (user equipment) in the cell performs the channel measurement through the pilot signals, which is used to calculate channel parameters such as CQI (channel status information) and to feed back the channel parameters to the base station. In the conventional MIMO technology, there are few transmitting antennas (generally no more than 8 antennas), thus each of the antennas transmits orthogonal pilot signals. However, in an MIMO system with large-scale antennas, with the increase of antennas, massive orthogonal pilot signals occupy a lot of time-frequency resources and reduce the actual system throughput. In view of this, the following method is provided in conventional technology.

The transmitting antennas are divided into multiple groups, the antennas in one transmit a same pilot signal, and the antennas in different groups transmit orthogonal pilot signals, i.e., the antennas in different groups transmit the pilot signals in a frequency-division mode. In the method, the same pilot signal is transmitted on the antennas in one group, which may save pilot overhead caused by multiple antennas, but the receiving end of the pilot signal cannot perform a complete channel measurement.

SUMMARY

In view of this, a channel measurement method for large-scale antennas and a user terminal thereof are provided in the present disclosure, to reduce pilot overhead and achieve integrity of channel measurement.

The technical solution is as follows.

In a first aspect, a channel measurement method for large-scale antennas is provided, where antennas at a base station side are divided into more than two groups in advance. The channel measurement method includes:

receiving, by a user terminal, pilot signals transmitted from different antennas in one at the base station side in a frequency-division mode and transmitted from antennas in different groups at the base station side in a time-division mode; and

performing a spatial channel estimation based on the pilot signals received from partial groups of antennas, to obtain a channel measurement result for all the antennas.

In a first possible implementation in the first aspect, the time-frequency resources each time occupied by each of the antennas in a sub-frame assigned to transmit the pilot signals are variable.

In a second possible implementation in the first aspect, performing the spatial channel estimation based on the pilot signals received from partial groups of antennas, to obtain the channel measurement result for all the antennas, includes:

wherein pilot signals y received by a mobile terminal from the partial groups of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas included in the partial groups of antennas;

based on an equation Ĥ=C′_(H)(C″_(H)+σ²diag(1/|x₁|², . . . , 1/|x_(K)|²))⁻¹y, obtaining a channel estimation value Ĥ for all the antennas, where Ĥ=(h₁, h₂, . . . , h_(M))^(T), M is the number of all the antennas, σ is noise power of a receiver in the user terminal, C′_(H) is column K in C_(H) corresponding to positions of the pilot signals received by the mobile terminal, C″_(H) cliff is row K and column K in C_(H) corresponding to positions of the pilot signals received by the mobile terminal, and the C_(H) is a channel correlation matrix.

Combined with the second possible implementation in the first aspect, in the case that the mobile terminal performs estimation of Ĥ for the first time, a preset initial value of the C_(H) is used; and

after Ĥ is obtained through each estimation, current C_(H) is updated based on a calculated Ĥ and is used for calculation to perform a next estimation of Ĥ, and the following updating equation is used to update the current C_(H) based on the calculated Ĥ: C _(H)=(1−a)C _(H) +aĤĤ ^(H), where Ĥ^(H) is conjugation of Ĥ, and α is a preset parameter.

In a third possible implementation in the first aspect, performing the spatial channel estimation based on the pilot signals received from the partial groups of antennas, to obtain the channel measurement result for respective antennas, includes:

wherein the pilot signals y received by a mobile terminal from the partial groups of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas included in the partial groups of antennas;

determining a channel value Ĥ′ for the K antennas in the partial groups of antennas, where

${{\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}};$

performing discrete Fourier transform on Ĥ′ at a point K, to obtain h₁;

padding (M−K) zero to an end of h₁ to obtain h₂ with M points; and

performing inverse discrete Fourier transform on h₂, to obtain a channel estimation value Ĥ for all the antennas.

In a fourth possible implementation in the first aspect, performing the spatial channel estimation based on the pilot signals received from the partial groups of antennas, to obtain the channel measurement result for respective antennas, includes:

where the pilot signals y received by a mobile terminal from the partial groups of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas included in the partial groups of antennas;

determining a channel value Ĥ′ for the K antennas in the partial groups of antennas, where

${{\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}};$

padding zero to positions in Ĥ′ to an antenna which does not transmit the pilot signals to obtain a matrix Ĥ″; and

performing spatial filtering on the Ĥ″ based on a preset spatial filter, to obtain a channel estimation value Ĥ for all the antennas.

In a second aspect, a user terminal is provided. The user terminal includes:

a pilot receiving unit, configured to receive pilot signals transmitted from different antennas in one group at a base station side in a frequency-division mode and transmitted from antennas in different groups at the base station side in a time-division mode; and

a channel estimation unit, configured to perform a spatial channel estimation based on the pilot signals received from partial groups of antennas, to obtain a channel measurement result of all the antennas.

In a first possible implementation in the second aspect, the pilot signals y received by a mobile terminal from the partial groups of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas included in the partial groups of antennas; and

the channel estimation unit includes:

a minimum mean square error estimation subunit, configured to, based on an equation Ĥ=C′_(H)(C″_(H)+σ²diag(1/|x₁|², . . . , 1/|x_(K)|²))⁻¹y, obtain a channel estimation value Ĥ for all the antennas, where Ĥ=(h₁, h₂, . . . , h_(M))^(T), M is the number of all the antennas, σ is noise power of a receiver in the user terminal, C′_(H) is column K in C_(H) corresponding to positions of the pilot signals received by the mobile terminal, C″_(H) is row K and column K in C_(H) corresponding to positions of the pilot signals received by the mobile terminal, and C_(H) is a channel correlation matrix.

Combined with the first possible implementation in the second aspect, the channel estimation unit further includes:

a correlation matrix determination subunit, configured to determine the C_(H) as a preset initial value in the case that the mobile terminal performs estimation for Ĥ for the first time; update a current C_(H) based on a calculated Ĥ after the minimum mean square error estimation subunit obtains Ĥ through estimation each time, to make the minimum mean square error estimation subunit use the current C_(H) to calculate for next estimation for Ĥ; where the following updating equation is used to update the current C_(H) based on the calculated Ĥ: C _(H)=(1−a)C _(H) +aĤĤ ^(H), where Ĥ^(H) is conjugation of Ĥ, and α is a preset parameter.

In a second possible implementation in the second aspect, the pilot signals y received by a mobile terminal from the partial groups of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas included in the partial groups of antennas; and

the channel estimation unit includes:

a channel value determination subunit, configured to determine a channel value Ĥ′ for the K antennas in the partial groups of antennas, where

${{\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}};$

a Fourier transform subunit, configured to perform discrete Fourier transform on Ĥ′ at a point K, to obtain h₁; and

an inverse Fourier transform subunit, configured to pad (M−K) zero to an end of the h₁ to obtain a matrix h₂ with M points, and perform inverse discrete Fourier transform on the h₂, to obtain a channel estimation value Ĥ for all the antennas.

In a third possible implementation in the second aspect, the pilot signals y received by a mobile terminal from the partial groups of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas included in the partial groups of antennas; and

the channel estimation unit includes:

a channel value determination subunit, configured to determine a channel value Ĥ′ for the K antennas in the partial groups of antennas, where

${{\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}};$

a zero padding subunit, configured to pad zero to positions in Ĥ′ to an antenna which does not transmit the pilot signals, to obtain a matrix Ĥ″; and

a spatial filtering subunit, configured to perform spatial filtering on Ĥ″ based on a preset spatial filter, to obtain a channel estimation value Ĥ for al the antennas.

It may be learned from above technical solution that, in the present disclosure, different antennas in one group at the base station side transmit pilot signals in the frequency-division mode, and antennas in different groups transmit the pilot signals in the time-division mode, which reduces pilot overhead caused by multiple antennas; the user terminal performs the spatial channel estimation based on the pilot signals received from the partial groups of antennas, to obtain the channel measurement result for all the antennas, thus a complete channel measurement is ensured.

In addition, the user terminal performs the spatial channel estimation without waiting for reception of pilot signals from all the antennas, to reduce delay for channel parameter feedback.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram of an MIMO system with large-scale antennas provided in the conventional technology;

FIG. 2 is a flow chart of a method according to an embodiment of the disclosure;

FIG. 3 is a diagram that each group antennas at a base station side transmit pilot signals according to an embodiment of the disclosure;

FIG. 4 is a diagram that one group of antennas transmit pilot signals according to the embodiment of the disclosure;

FIG. 5 is a diagram of a first structure for a user terminal according to an embodiment of the disclosure;

FIG. 6 is a diagram of a second structure for the user terminal according to an embodiment of the disclosure;

FIG. 7 is a diagram of a third structure for the user terminal according to an embodiment of the disclosure; and

FIG. 8 is a graph that MSE varies with SNR under different values of K according to an embodiment of the disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

To make the objective, technical solution and merits of the disclosure more clearly, the disclosure is described in detail in conjunction with drawings and specific embodiments hereinafter.

The method according to the embodiments of the disclosure is shown in FIG. 2. The method mainly includes steps at two sides, i.e., steps at a base station side and steps at a terminal side.

In step 201, antennas at the base station side are divided into more than two groups in advance, the different antennas in one group transmit pilot signals in a frequency-division mode, and the antennas in different groups transmit the pilot signals in a time-division mode.

In step 202, a UE side performs a spatial channel estimation based on the pilot signals received from partial groups of the antennas, to obtain a channel measurement result for each of the antennas.

Each of the steps is described in detail in conjunction with the embodiments hereinafter. Supposing that there are M antennas at the base station side, and the antennas at the base station side are divided into N groups, where N is an integer greater than 1, each of the groups includes the K antennas, and K=M/N. To make spatial pilots even and a channel estimation error small, preferably, the antennas in groups are distributed at a regular interval.

In the step 201, the antennas in different groups transmit the pilot signals in the time-division mode. It may be set that the antennas in one group transmit the pilot signals every n sub-frames, where n is an integer greater than or equal to N, thus antennas in different groups transmit the pilot signals every n sub-frames, and pilot signals from different antennas are transmitted in different sub-frames. In the group, different antennas occupy different time-frequency resources. In addition, preferably, each of the antennas may use different time-frequency resources to transmit a pilot signal in the adopted sub-frames every time. Supposing that group 1 of the antennas transmit pilot signals at the first sub-frame, the (n+1)th sub-frame, the (2n+1)th sub-frame . . . , antenna 11 in the group 1 of antennas respectively transmits the pilot signals in the first sub-frame by using the first time-frequency resource, in the (n+1)th sub-frame by using the second time-frequency resource, and in the (2n+1)th sub-frame by using the third time-frequency resource . . . as long as there is an appointment which ensures that the time-frequency resources occupied by the antennas in one group do not conflict with each other.

For example, supposing there are eighteen antennas, which are divided into three groups with six antennas each group. Each group of antennas transmit pilot signals every n sub-frames, the antennas in group 1, the antennas in group 2 and the antennas in group 3 transmit the pilot signals as shown in FIG. 3, and boxes with bold lines represent the sub-frames where the groups of antennas transmit the pilot signals.

The antennas in the group 1 are taken as an example. The six antennas in the same group respectively transmit the pilot signals by using the time-frequency resources shown as 6 small squares in the sub-frame. As shown in FIG. 4, the solid squares in the figure represent the time-frequency resources occupied by the antenna 11, and the antenna 11 may occupy different time-frequency resources in different sub-frames. For example, in the first sub-frame, the antenna 11 occupies the time-frequency resource in the first row and first column; in the (n+1)th sub-frame, the antenna 11 occupies the time-frequency resource in the second row and first column. That is, the time-frequency resources occupied by each of the antennas in the sub-frame assigned to transmit the pilot signals are variable each time.

In the step 202, the UE side may perform the spatial channel estimation based on the pilot signals received from partial groups of antennas, to obtain the channel measurement result for each of the antennas. A case that the UE side uses the pilot signals received from the antennas in one group is taken as example, for description of the spatial channel estimation. Supposing that, in a sub-frame, K orthogonal pilots transmitted by the antennas in the group 1 are x=(x₁, x₂, . . . , x_(K))^(T), pilot signals received by the UE from the antennas in the group 1 are y=(y₁, y₂, . . . , y_(K))^(T), and the spatial channel estimation performed in the disclosure is to obtain a channel measurement result Ĥ, where Ĥ is a channel estimation value for M antennas, Ĥ=(h₁, h₂, . . . , h_(K))^(T), and M may be the quantity of the antennas. The following three ways may be used for the spatial channel estimation according to the embodiments of the disclosure, but the disclosure is not limited hereto.

A first way is a spatial LMMSE (linear minimum mean square error) estimation.

The UE obtains Ĥ based on the received pilot signals y, pilot signals x transmitted by the antennas in the group 1 and a channel correlation matrix C_(H), i.e., the UE estimates channel value of each of the antenna ports based on a correlation among the antennas. In this way, Ĥ may be calculated based on the following equation (1): Ĥ=C′ _(H)(C″ _(H)+σ²diag(1/|x ₁|², . . . ,1/|x _(K)|²))⁻¹ y  (1) Where σ is noise power of a receiver at the UE side,

${C_{H}^{\prime} = \begin{pmatrix} c_{11} & c_{1,{({K + 1})}} & L & c_{1,M} \\ c_{21} & O & \; & \; \\ M & \; & O & \; \\ M & \; & \; & \; \\ c_{M,1} & \; & \; & c_{M,M} \end{pmatrix}_{M \times K}},$ i.e., in C_(H), column K corresponding to positions of the pilot signals received by the UE is taken.

${C_{H}^{''} = \begin{pmatrix} c_{11} & c_{1,{({K + 1})}} & L & c_{1,M} \\ c_{{({K + 1})},1} & O & \; & \; \\ M & \; & O & \; \\ c_{M,1} & \; & \; & c_{M,M} \end{pmatrix}_{K \times K}},$ in C_(H), row K and column K corresponding to the positions of the pilot signals received by the UE is taken.

When the UE performs estimation of Ĥ for the first time, a preset initial value of C_(H) is used. After Ĥ is obtained through each calculation, the current C_(H) is updated based on the calculated Ĥ and is used for calculation to perform a next estimation of Ĥ. That is, an updated C_(H) is used for calculation each time. The following updating equation (2) is used to update the current C_(H) based on Ĥ:

$\begin{matrix} {C_{H} = {{{\left( {1 - \alpha} \right)C_{H}} + {\alpha\begin{pmatrix} O & \; & O \\ \; & {{\hat{h}}_{i}{\hat{h}}_{j}^{*}} & \; \\ O & \; & O \end{pmatrix}}} = {{\left( {1 - \alpha} \right)C_{H}} + {\alpha\;\hat{H}{\hat{H}}^{H}}}}} & (2) \end{matrix}$

where Ĥ^(H) is conjugation of Ĥ, α is a preset parameter value, which may be set based on the demands for update accuracy and speed of the correlation matrix, and is generally set at the order of magnitude of 0.01.

A second way is a spatial DFT (discrete Fourier transform) interpolation algorithm.

At first, the channel value Ĥ′ of the K antennas transmitting the pilot signals is determined based on the following equation:

${\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}$

DFT is performed on Ĥ′ at K points, to obtain h₁, and h₁=DFT {Ĥ′}.

Then, (M−K) zero is padded to an end of h₁ to obtain h₂, where the length of h₂ is M points, and h₂=(h₁, 0, . . . , 0).

At last, IDFT (inverse discrete Fourier transform) is performed on h₂, and Ĥ is obtained, and Ĥ=IDFT{h₂}.

A third way is spatial filtering.

At first, the channel value Ĥ′ of the K antennas transmitting the pilot signals is determined based on the following equation:

${\hat{H}}^{\prime} = {\left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}.}$

Zero is padded to positions in Ĥ′ corresponding to the antennas which do not transmit the pilot signals to obtain Ĥ″:

${\hat{H}}^{''} = {\left( {\frac{y_{1}}{x_{1}},0,{\ldots\mspace{14mu}\frac{y_{2}}{x_{2}}},0,\ldots\mspace{14mu},\frac{y_{K}}{x_{K}},0,\ldots}\mspace{14mu} \right)^{T}.}$

Spatial filtering is performed on Ĥ″ based on a preset spatial filter, to obtain Ĥ. The spatial filter may adopt methods such as conventional window function, and has a cut-off frequency set as

${\frac{f_{s}}{2}\frac{K}{M}},$ where f_(s) is a sampling frequency of the spatial filter, i.e., a normalized cut-off frequency of the spatial filter is

$\frac{K}{M}.$

According to the embodiment of the disclosure, the mobile terminal performs the channel estimation without waiting for reception of pilot signals from all the antennas, and may obtain the channel measurement result for the respective antennas through spatial channel estimation after receiving the pilot signals from partial groups of antennas. For example, the channel measurement result for the respective antennas may be obtained based on the pilot signals from the group of antennas after the pilot signals from one group of antennas are received. After the channel estimation result is obtained, channel parameters may be calculated and fed back, which greatly reduces feedback delay.

The above is description of the method according to the embodiment of the disclosure, and the following is description of a user terminal according to the embodiment of the disclosure.

FIG. 5 is a diagram of structure for the user terminal according to the embodiment of the disclosure. As shown in FIG. 5, the user terminal includes a pilot receiving unit 00 and a channel estimation unit 10.

The pilot receiving unit 00 is used for receiving pilot signals transmitted by a base station side. Specifically, the base station side transmits the pilot signals in this way: antennas at the base station side are divided into more than two groups in advance, the antennas in one group transmit the pilot signals in a frequency-division mode, i.e., transmit orthogonal pilots, and the antennas in different groups transmit the pilot signals in a time-division mode.

The channel estimation unit 10 performs a spatial channel estimation based on pilot signals received from partial groups of the antennas, to obtain a channel measurement result for each of the antennas. According to the embodiment of the disclosure, the mobile terminal performs the channel estimation without waiting for reception of pilot signals from the respective antennas, and may obtain the channel measurement result for the respective antennas through spatial channel estimation after receiving the pilot signals from partial groups of antennas.

According to the embodiment of the disclosure, the channel estimation unit 10 may use multiple methods for spatial channel estimation. Corresponding to different methods, the channel estimation unit 10 may be implemented in different ways, which include but are not limited to the following.

A first implementation is shown in FIG. 5. The channel estimation unit 10 includes a minimum mean square error estimation subunit 11 which uses a spatial LMMSE (linear minimum mean square error) estimation.

It is supposed that pilot signals y received by the mobile terminal from the partial groups (supposing that there is one group of antennas, and K is the number of the antennas in the group) of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas in the partial groups.

Based on an equation Ĥ=C′_(H)(C″_(H)+σ²diag(1/|x₁|², . . . , 1/|x_(K)|²))⁻¹y, the minimum mean square error estimation subunit 11 may obtain channel estimation value Ĥ for respective channels, Ĥ=(h₁, h₂, . . . , h_(M))^(T), where M is the number of the all antennas, σ is noise power of a receiver at a UE side, C′_(H), is column K in C_(H) corresponding to positions of the pilot signals received by the mobile terminal, C″_(H) is row K and column K in C_(H) corresponding to positions of the pilot signals received by the mobile terminal, and C_(H) is a channel correlation matrix.

Since the minimum mean square error estimation subunit 11 uses the channel correlation matrix C_(H) when performing estimation for Ĥ, there is a preferred implementation, i.e., C_(H) may be updated iteratively based on each estimated Ĥ. In this case, the channel estimation unit may further include a correlation matrix determination subunit 12.

The correlation matrix determination subunit 12 is configured to determine C_(H) as a preset initial value in the case that the mobile terminal performs estimation for Ĥ for the first time; update a current C_(H) based on a calculated Ĥ if after the minimum mean square error estimation subunit 11 obtains Ĥ through calculation each time, to make the minimum mean square error estimation subunit 11 use the current C_(H) to calculate for next estimation for Ĥ The following equation is used when the current C_(H) is updated through the calculated Ĥ: C _(H)=(1−a)C _(H) +aĤĤ ^(H), where Ĥ^(H) is conjugation of Ĥ, and α is a preset parameter.

A second implementation is shown in FIG. 6. The channel estimation unit 10 uses a spatial DFT (discrete Fourier transform) interpolation algorithm. The channel estimation unit 10 includes a channel value determination subunit 21, a Fourier transform subunit 22, a zero padding processing subunit 23 and an inverse Fourier transform subunit 24.

It is supposed that pilot signals y received by the mobile terminal from the partial groups of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas included in the partial groups.

At first, the channel value determination subunit 21 determines channel value Ĥ′ for the K antennas in the partial groups of antennas,

${\hat{H}}^{\prime} = {\left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}.}$ The Fourier transform subunit 22 performs DFT on Ĥ′ at a point K, to obtain h₁. Then, the zero padding processing subunit 23 pads (M−K) zero to the end of the h₁ to obtain a matrix h₂ with M points. At last, the inverse Fourier transform subunit 24 performs inverse discrete Fourier transform on the h₂, and channel estimation value Ĥ for all the antennas are obtained.

A third implementation is shown in FIG. 7. The channel estimation unit 10 uses spatial filtering. The channel estimation unit 10 includes a channel value determination subunit 31, a zero padding subunit 32 and a spatial filtering subunit 33.

It is supposed that pilot signals y received by the mobile terminal from the partial groups of antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), where K is the number of the antennas included in the partial groups.

At first, the channel value determination subunit 31 determines channel value Ĥ′ for the K antennas in the partial groups of antennas,

${\hat{H}}^{\prime} = {\left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}.}$ Then, the zero padding subunit 32 pads zero to positions in Ĥ′ corresponding to antennas which do not transmit the pilot signals, to obtain Ĥ″. At last, the spatial filtering subunit 33 performs spatial filtering on Ĥ″ based on a preset spatial filter, to obtain channel estimation value Ĥ for all the antennas. The spatial filter may adopt methods such as conventional window function, and has a cut-off frequency set as

${\frac{f_{s}}{2}\frac{K}{M}},$ where f_(s) is a sampling frequency of the spatial filter, i.e., a normalized cut-off frequency of the spatial filter is

$\frac{K}{M}.$

According to the method and user terminal provided in the disclosure, pilot overhead caused by multiple antennas may be effectively reduced, and errors in channel estimation may be effectively controlled, which is verified based on a set of experimental data.

Supposing that the number M of all antennas is 600, experiments are conducted in cases that the numbers K of received orthogonal pilots are 600, 120, 60, 30, 20 and 12 respectively. Supposing that an interval in time-domain is 5 ms, proportions of time-frequency resources occupied by the orthogonal pilots are shown in Table 1.

TABLE 1 K 600 120 60 30 20 12 intervals for 1 5 10 20 30 50 respective antennas to transmit orthogonal pilots proportions of 71.4% 14.3% 7.14% 3.6% 2.4% 1.4% time-frequency resources occupied by the orthogonal pilots

It may be seen from Table 1 that, in a case that all the 600 antennas transmit orthogonal pilots, occupied time-frequency resources exceed 70%, and remaining time-frequency resources available for data transmission are less than 30%. In a case that the orthogonal pilots are transmitted by means of grouping and round-robin according to the embodiments of the disclosure, supposing that the orthogonal pilots are transmitted 60 antennas, the time-frequency resources occupied by the orthogonal pilots are only about 7%, and time-frequency resources available for data transmission are greatly improved.

FIG. 8 is a graph that a channel estimation error (MSE) varies with SNR (signal noise ratio) under different values of K according to an embodiment of the disclosure. It may be seen that, even if the value of K is 60, the channel estimation error is below 0.1, and the accuracy is relatively high. It may be also seen from the figure that, the more groups there are, i.e., the smaller the value of K is, and the larger the channel estimation error, thus it is required to make a tradeoff between consumption of time-frequency resources and the channel estimation error to elect the value of K.

According to the embodiments of the disclosure, it should be understood that, the described system, device and methods may be implemented in other ways. For example, the described device embodiment is merely exemplary. The division of the units is merely based on logical functions, and the units may be divided with other approaches in practice.

In addition, individual functional units according to the embodiments of the disclosure may be integrated in one processing unit, or the units may exist separately, or two or more units may be integrated in one unit. The foregoing integrated units may be realized in a form of hardware, or realized in a form of functional units with hardware and software.

The integrated unit realized in the form of software function unit may be stored in a computer readable storage medium. The software function unit is stored in a storage medium and includes some instructions for instructing a computer device (may be personal computer, server, network device, etc.) or a processor (processor) to implement part of the steps of the method according to the embodiments of the disclosure. The foregoing storage medium includes various media that can store program codes, for example, USB disk, mobile hard disk drive, read-only memory (Read-Only Memory, ROM), random access memory (Random Access Memory, RAM), magnetic disk, or optical disk and the like.

Only preferred embodiments of the disclosure are described above, which are not to limit the scope of the disclosure. Any change, substitution or modification within the spirit and principle of the disclosure fall within the protection scope of the application. 

What is claimed is:
 1. A channel measurement method for large-scale antennas, wherein antennas at a base station side are divided into more than two groups and the channel measurement method comprises: receiving, by a user terminal, pilot signals transmitted from the antennas at the base station side, wherein the antennas in one group transmit the pilot signals in a frequency-division mode, the antennas in different groups transmit the pilot signals in a time-division mode, and the antennas in different groups transmit the pilot signals in different sub-frames; and performing a spatial channel estimation based on the pilot signals received from partial groups of the antennas, to obtain a channel measurement result for all the antennas.
 2. The channel measurement method according to claim 1, wherein time-frequency resources each time occupied by each of the antennas in a sub-frame assigned to transmit the pilot signals are variable.
 3. The channel measurement method according to claim 1, wherein performing the spatial channel estimation based on the pilot signals received from the partial groups of the antennas, to obtain the channel measurement result for all the antennas, comprises: wherein the pilot signals Y received by the user terminal from the partial groups of the antennas is y=(y₁, y₂, . . . , y_(K))^(T), and the pilot signals x transmitted by the partial groups of the antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), wherein K is the number of the antennas comprised in the partial groups of antennas; and based on an equation Ĥ=C′_(H)(C″_(H)+σ²diag(1/|x|², . . . , 1/|x_(K)|²))⁻¹y, obtaining a channel estimation value Ĥ for all the antennas, wherein H=(h₁, h₂, . . . , h_(M))^(T), M is the number of all the antennas, σ is noise power of a receiver in the user terminal, C′_(H) is column K in C_(H) corresponding to positions of the pilot signals received by the user terminal, C″_(H) is row K and column K in C_(H) corresponding to positions of the pilot signals received by the user terminal, and the C_(H) is a channel correlation matrix.
 4. The method according to claim 3, wherein in the case that the user terminal performs estimation of Ĥ for the first time, a preset initial value of the C_(H) is used; and after Ĥ is obtained through each estimation, current C_(H) is updated based on a calculated Ĥ and is used for calculation to perform a next estimation of Ĥ, and the following updating equation is used to update the current C_(H) based on the calculated Ĥ: C_(H)=(1−a)C_(H)+aĤĤ^(H), wherein Ĥ^(H) is a conjugate transpose of Ĥ, and α is a preset parameter.
 5. The method according to claim 1, wherein performing the spatial channel estimation based on the pilot signals received from the partial groups of the antennas, to obtain the channel measurement result for all the antennas, comprises: the pilot signals y received by the user terminal from the partial groups of the antennas is y=(y₁, y₂, . . . , y_(K))^(T), and the pilot signals x transmitted by the partial groups of the antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), wherein K is the number of the antennas comprised in the partial groups of the antennas; determining channel values Ĥ′ for the K antennas in the partial groups of the antennas, wherein ${{\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}};$ performing discrete Fourier transform on Ĥ′ at K points, to obtain h₁; padding (M−K) zero to an end of h₁ to obtain h₂ with M points; and performing inverse discrete Fourier transform on h₂, to obtain a channel estimation value Ĥ for all the antennas.
 6. The method according to claim 1, wherein performing the spatial channel estimation based on the pilot signals received from partial groups of the antennas, to obtain the channel measurement result for all the antennas, comprises: wherein the pilot signals y received by the user terminal from the partial groups of the antennas is y=(y₁, y₂, . . . , y_(K))^(T), and pilot signals x transmitted by the partial groups of the antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), wherein K is the number of the antennas comprised in the partial groups of the antennas; determining channel values Ĥ′ for the K antennas in the partial groups of the antennas, wherein ${{\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}};$ padding zero to positions in Ĥ′ corresponding to antennas which do not transmit the pilot signals to obtain a matrix Ĥ″; and performing spatial filtering on the Ĥ″ based on a preset spatial filter, to obtain a channel estimation value Ĥ for all the antennas.
 7. A user terminal, comprising a receiver and a processor, wherein: the receiver receives pilot signals transmitted from antennas at a base station side, wherein the antennas in one group transmit the pilot signals in a frequency-division mode, the antennas in different groups transmit the pilot signals in a time-division mode, and the antennas in different groups transmit the pilot signals in different sub-frames; and the processor performs a spatial channel estimation based on the pilot signals received from partial groups of the antennas, to obtain a channel measurement result of all the antennas.
 8. The user terminal according to claim 7, wherein the pilot signals y received by the user terminal from the partial groups of the antennas is y=(y₁, y₂, . . . , y_(K))^(T), and the pilot signals x transmitted by the partial groups of the antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), wherein K is the number of the antennas comprised in the partial groups of the antennas; and the processor, based on an equation Ĥ=C′_(H) (C″_(H)+σ²diag(1/|x₁|², . . . , 1/|x_(K)|²))⁻¹y, obtains channel estimation values Ĥ for all the channels, wherein Ĥ=(h₁, h₂, . . . , h_(M))^(T), M is the number of all the antennas, σ is noise power of a receiver in the user terminal, C′_(H) is column K in C_(H) corresponding to positions of the pilot signals received by the user terminal, C″_(H) is row K and column K in C_(H) corresponding to positions of the pilot signals received by the user terminal, and C_(H) is a channel correlation matrix.
 9. The user terminal according to claim 8, wherein the processor determines the C_(H) as a preset initial value in the case that the user terminal performs estimation for Ĥ for the first time; and updates a current C_(H) based on a calculated Ĥ after the minimum mean square error estimation subunit obtains Ĥ through calculation each time, to make the minimum mean square error estimation subunit use the current C_(H) to calculate for next estimation for Ĥ; wherein the following updating equation is used to update the current C_(H) based on the calculated Ĥ: C_(H)=(1−a)C_(H)+aĤ Ĥ^(H), wherein Ĥ^(H) is a conjugate transpose of Ĥ, and α is a preset parameter.
 10. The user terminal according to claim 7, wherein the pilot signals y received by the user terminal from the partial groups of the antennas is y=(y₁, y₂, . . . , y_(K))^(T), and the pilot signals x transmitted by the partial groups of the antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), wherein K is the number of the antennas comprised in the partial groups of the antennas; and the processor determines channel values Ĥ′ for the K antennas in the partial groups of the antennas, ${{\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}};$ performs discrete Fourier transform on Ĥ′ at K points, to obtain h₁; and pads (M−K) zero to an end of the h₁ to obtain a matrix h₂ with M points, and perform inverse discrete Fourier transform on the h₂, to obtain a channel estimation value Ĥ for all the antennas.
 11. The user terminal according to claim 7, wherein the pilot signals y received by the user terminal from the partial groups of the antennas is y=(y₁, y₂, . . . , y_(K))^(T), and the pilot signals x transmitted by the partial groups of the antennas at the base station side is x=(x₁, x₂, . . . , x_(K))^(T), wherein K is the number of the antennas comprised in the partial groups of the antennas; and the processor determines a channel value Ĥ′ for the K antennas in the partial groups of the antennas, ${{\hat{H}}^{\prime} = \left( {\frac{y_{1}}{x_{1}},\frac{y_{2}}{x_{2}},\ldots\mspace{14mu},\frac{y_{K}}{x_{K}}} \right)^{T}};$ pads zero to positions in Ĥ′ corresponding to antennas which do not transmit the pilot signals, to obtain a matrix Ĥ″; and performs spatial filtering on Ĥ″ based on a preset spatial filter, to obtain a channel estimation value Ĥ for all the antennas. 